A second order virtual node method for elliptic problems with interfaces and irregular domains in three dimensions
نویسندگان
چکیده
We present a numerical method for the variable coefficient Poisson equation in three dimensional irregular domains and with interfacial discontinuities. The discretization embeds the domain and interface into a uniform Cartesian grid augmented with virtual degrees of freedom to provide accurate treatment of jump and boundary conditions. The matrix associated with the discretization is symmetric positive definite and equal to the standard 7-point finite difference stencil away from embedded interfaces and boundaries. Numerical evidence suggests second order accuracy in the L∞-norm. Our approach improves the treatment of Dirichlet and jump constraints in the recent work of Bedrossian et al. [11] and provides novel aspects necessary for problems in three dimensions. Specifically, we construct new constraint-based Lagrange multiplier spaces that significantly improve the conditioning of the associated linear system of equations; we provide a method for sub-cell polyhedral approximation to the zero isocontour surface of a level set needed for three dimensional embedding; and we show that the new Lagrange multiplier spaces naturally lead to a class of easy-to-implement multigrid methods that achieve near-optimal efficiency, as shown by numerical examples. For the specific case of a continuous Poisson coefficient in interface problems, we provide an expansive treatment of the construction of a particular solution that satisfies the value jump and flux jump constraints. As in [11], this is used in a discontinuity removal technique that yields the standard 7-point Poisson stencil across the interface and only requires a modification to the right-hand side of the linear system.
منابع مشابه
A second order virtual node method for elliptic problems with interfaces and irregular domains
We present a second order accurate, geometrically flexible and easy to implement method for solving the variable coefficient Poisson equation with interfacial discontinuities or on irregular domains, handling both cases with the same approach. We discretize the equations using an embedded approach on a uniform Cartesian grid employing virtual nodes at interfaces and boundaries. A variational me...
متن کاملVirtual Node Algorithms for Stokes Interface Problems
of the Dissertation Virtual Node Algorithms for Stokes Interface Problems by Diego C. Assêncio Doctor of Philosophy in Physics University of California, Los Angeles, 2012 Professor Giovanni Zocchi, Chair We present two numerical methods for the solution of the Stokes equations designed to handle both interfacial discontinuities, geometrically irregular flow domains and discontinuous fluid prope...
متن کاملMatched interface and boundary (MIB) method for elliptic problems with sharp-edged interfaces
Elliptic problems with sharp-edged interfaces, thin-layered interfaces and interfaces that intersect with geometric boundary, are notoriously challenging to existing numerical methods, particularly when the solution is highly oscillatory. This work generalizes the matched interface and boundary (MIB) method previously designed for solving elliptic problems with curved interfaces to the aforemen...
متن کاملOn the Solution of the Elliptic Interface Problems by Difference Potentials Method
Designing numerical methods with high-order accuracy for problems in irregular domains and/or with interfaces is crucial for the accurate solution of many problems with physical and biological applications. The major challenge here is to design an efficient and accurate numerical method that can capture certain properties of analytical solutions in different domains/subdomains while handling ar...
متن کاملHigh order matched interface and boundary method for elliptic equations with discontinuous coefficients and singular sources
This paper introduces a novel high order interface scheme, the matched interface and boundary (MIB) method, for solving elliptic equations with discontinuous coefficients and singular sources on Cartesian grids. By appropriate use of auxiliary line and/or fictitious points, physical jump conditions are enforced at the interface. Unlike other existing interface schemes, the proposed method disas...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- J. Comput. Physics
دوره 231 شماره
صفحات -
تاریخ انتشار 2012